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Abstract 

We calculate couplings of arbitrary order from correlation functions among 
twisted strings, using conformal field theory. Twisted strings arise in heterotic string 
compactified on orbifolds yielding matter fields in the low energy limit. We calculate 
completely the classical and the quantum amplitude including normalization, up to 
a contribution from Kahler potential. The classical action has saddle points which 
are interpreted as worldsheet instantons described by metastable untwisted strings, 
formed by twisted strings distributed at certain fixed points. This understanding 
generalizes the area rule, in the case that the locations of twisted strings do not 
form a polygon, and provides a general rule for calculating these kinds of instanton 
corrections. An interpretation of couplings involving linearly combined states is 
given, which commonly appear in non-prime order orbifolds. The quantum part of 
the amplitude is given by ratios of gamma functions with order one arguments. 
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1 Introduction 



Superstring theory is a promising candidate for unified theory including gravity. Heterotic 
orbifold construction is one of interesting constructions for four dimensional string models 
[HE]. (See also for resent works Ref. [31 H] and for review [5].) One can solve equation of 
motion of string on the orbifold background, and geometrical picture is clear in heterotic 
orbifold models. Thus, several aspects can be computed and can be understood from the 
geometrical viewpoint. 

Heterotic orbifold models have modes localized at fixed points, that is, twisted strings. 
3-point couplings as well as 4-point couplings of these localized modes have been computed 
analytically [6l O |8l El |10], and the size of Yukawa coupling Y is obtained as F ~ e~^, 
where A denotes naively the area of triangle corresponding to three fixed points of twisted 
strings. This will be clarified more and generalized in this paper. This aspect is quite 
interesting from the phenomenological viewpoint. One can obtain suppressed Yukawa 
couplings when twisted strings are localized far away from each other. That is, one could 
explain the hierarchy of quark and lepton masses as well as their mixing angles when they 
are localized at different places. 

We have to study selection rules for allowed couplings in order to examine whether 
realistic fermion masses and mixing angels can be realized from string theory. The space 
group selection rule [H [H] constrain allowed Yukawa couplings rather strongly. For ex- 
ample, on prime order orbifolds off-diagonal Yukawa couplings are not allowed, and we 
can not obtain realistic mixing angles by using only 3-point couplings with the minimum 
number of Higgs fieldsQ On non-prime order orbifolds, off-diagonal Yukawa couplings are 
allowed pT| and possibilities for leading to realistic quark and lepton masses and mixing 
angles have been studied \13\. However, to realize fermion masses and mixing angles in 
string theory is still a challenging issue. 

In this paper, we study generic higher order couplings than renormalizable couplings 
in heterotic orbifold models. Higher dimensional operators become effective Yukawa cou- 
plings after symmetry breaking. Suppose that we have a coupling of type FfH Yli 4>i 
in the superpotential of effective field theory, where F and / are chiral matter fields 
corresponding quarks and leptons, H denotes electroweak Higgs superfields and 0, corre- 
spond to several heavy modes. When all scalar components of the superfields (pi develop 
their vacuum expectation values (VEVs), this higher dimensional operator becomes a 
Yukawa coupling among chiral fermions F and / and the electroweak Higgs fields H. 
Thus, there is a possibility for deriving quark/lepton masses and mixing angles through 
this type of symmetry breaking, but by use of not only 3-point couplings. Indeed, such 
possibility has been examined in explicit models [2],[3lll]. Therefore, it is important to 
study selection rules of allowed higher order couplings and compute magnitude of allowed 

^In non-factorizable orbifold models, off-diagonal Yukawa couplings are allowed, but it is still difficult 
to derive realistic Yukawa matrices [T^] . 
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couplings^ When 0j correspond to locahzed modes on orbifold fixed points, the above ef- 
fective Yukawa coupling may correspond to a 3-point coupling on a Calabi-Yau manifold, 
where orbifold singularities are smoothed by the VEVs of 0^. Thus, calculations of higher 
order couplings on the orbifold are also important from the viewpoint of calculations of 
3-point couplings on the Calabi-Yau manifold around the orbifold limit. 

We compute magnitudes of L-point couplings. Generic aspects of L-point couplings 
heterotic orbifold models have been obtained in [8] . Here we apply it to concrete heterotic 
orbifold models. Similar calculation has been carried out for generic L-point couplings in 
intersecting D-brane models [T7]FI 

The paper is organized as follows. In section 2, we give a brief review on heterotic 
orbifold models in order to fix our notation. Then, we study the selection rule due to 
discrete i?-symmetry and the space group. In section 3, we compute classical contributions 
of L-point couplings. Their quantum parts are calculated in section 4. In section 5, we 
consider normalization of correlation functions. Section 6 is devoted to conclusion and 
discussion. In appendix, we give useful formulae for hypergeometric functions and their 
multivariable generalizations. 



2 Setup 

2.1 Twisted strings and their vertex operators 

The heterotic string theory consists of lOD right-moving superstring and 26D left-moving 
bosonic string. For the common ten (bosonic) dimensions, we consider the background 
with our 4D space-time and 6D orbifold. The other 16D left-moving bosonic string cor- 
respond to a gauge part. A 6D orbifold is a division of a 6D torus by a twist 6*, 
while is obtained as R^/A, where A is a 6D lattice. The twist 6 must be an auto- 
morphism of the lattice A, and its eigenvalues are diag(e^'^*''^ , e^'^*''^, e^'^^^^) in the complex 
basis Zi {i = 1,2,3). We mainly concentrate ourselves to the case that is factorizable 
as T'^ X T"^ X T^. To preserve 4D N = 1 supersymmetry (SUSY), they must satisfy the 
following condition, 

'ni+'n2 + V3= integer, (1) 

where r]i is not integer for each i = 1,2,3. 

The twisted string is a closed string up to orbifold identification 

Z{e^^'z,e-^^'z) =e''Z{z,z)+v, v e A, (2) 

where A is the above lattice (in the complex basis) defining the orbifold. It makes sense 
to restrict the phase to be —1 < kr]i < 1. Its zero mode satisfies the same condition, 

^ It would also be useful to study non-Abelian flavor symmetries [Mj [T^ and accidental global sym- 
metries |16j in string models. 

■^See for 3-point couplings in intersecting D-brane models \T8 \ [T9 l [20 l I21j. 
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and it is called a fixed point on the orbifold. The fixed point can be represented by the 
corresponding space group element, {9^,v). Note that the fixed point {O^^v) is equivalent 
to {e^,v + (1 - ^'=)A). They belong to the same conjugacy class. The sector with /c = 
corresponds to the so-called untwisted sector. 

The local operator called the twist operator ak{z) takes into account the nontrivial 
boundary condition ([2]) by inducing a branch point at z with the order k/N on the world- 
sheet, but the theory remains local by moding out by orbifold projection. The ground 
state corresponding to the twisted string on the fixed point {9^,v) is generated from the 
untwisted ground state |0) by the twist field \ak) = o"fc(0,0)|0). These twist fields have 
the operator product expansions (OPEs) 

9Z(^)orfc(0,0)~^'=/^-Vfc(0,0), (3) 
dZ{z)a,{0,0)r.z-'/''fUO,0), (4) 

which are understood as the most singular parts in the mode expansion for < ^ < 1. 
For the other region — 1 < < 0, we have the corresponding relations by replacing k 
with N — k. Also we have similar expressions for Z. Their conformal weights for the 
holomorphic and anti-holomorphic parts are 

_ _1A;/A;\ 

thus inducing a shift of zero point energy. 

Each 6''^-twisted sector has several ground states, that is, twist fields corresponding 
to several fixed points under 9^ twist. When we specify the fixed point /, we denote 
CT/^fe. Also we use the notation where {9^,v) denotes the space group element 

corresponding to the fixed point / under 9'^ twist. 

On non-prime order orbifolds, fixed points under higher twist 6**^ {k > 1) are not 
always fixed under 9 or twist fields (T(^gk^y^ are not always eigenstates of the twist 9. To 
make eigenstates, we have to take the following linear combinations [22| [TT]. 



where 7 = e^'^*^/'" with integer i to be determined by gauge quantum numbers and internal 
momenta. This linear combination may include twist fields corresponding to fixed points, 
which belong to the same conjugacy class. 

We consider the covariant quantization with the explicit conformal and superconformal 
ghosts. It is convenient to bosonize right-moving fermionic string and write bosonized 
degrees of freedom by H^{z). In the bosonized formulation, untwisted massless modes 
have momenta pt for ^ = ' ' ' 5 5, which are quantized on the 5*0(10) weight lattice. The 
space-time boson and fermion correspond to SO (10) vector and spinor, respectively. The 
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compact space corresponds to 5*0(6). The twisted sector has shifted SO (6) momenta, 
T^i = Pi + krji, which are often called if- moment a. 

A bosonic massless state has the corresponding vertex operator, 



J^iJnHtJP'x'ikXh) 



(7) 



i=l 



naturally in the (— l)-picture, where is the bosonized ghost, P^X^ corresponds to the 
gauge part and kX corresponds to 4D part. Here, d"^'Zi and d^^Zi denote oscillators 
for the left-mover, and Mi and Hi are oscillator numbers, which these massless modes 
include. Similarly, we can write massless modes corresponding to space-time fermions as 



F 1 = e-5<^]^(5'"'Z,)^"(9'^-Z0^'e^^''^^*e^^'^'e^'=^aJ],^. 



in the (— |)-picture. We understand that the H fields contains the four dimensional 

spin field. The if-momenta for space-time fermion and boson, r\ and in the same 
supersymmetric multiplet are related each other 

r, = rK^ + (l,l,l,l,l)/2, (9) 

that is, (1, 1, 1, 1, l)/2 corresponds to the ii-momentum of unbroken 4D space-time SUSY 
charge. To each vertex operator, we have to include overall normalization 

M r , . 1 /o n M 

(10) 



1/2 



(m, - 1)! 



1/2 



mi 



1)! 



from the state-operator mapping or the unitarity relation [23j . The closed string coupling 
is expressed in terms of ten dimensional gauge and gravitational couplings 



9c 



477 



2^' 



'111 



Thus, including one more field suppresses the corresponding coupling by one inverse mass 
dimension (9(a'^''^) as we expect. We have omitted the two-cocycles, which determine the 
overall sign [24j. 

We calculate a correlation function among L twisted matter fields including two space- 
time fermions on the Zjy orbifold, along the lines [3, [3 [9l [17]. It yields a higher order 
coupling in the zero momentum limit fcj 0. Since the background has the superconfor- 
mal ghost charge 2, the correlation function is of the form. 



„ L L 
lim e"^ / : ccV-x{,Zx,Zi) :: ccK_i(^2,^2) ccV_i(z3,53) : J]^ : Vq^Zi^z^ : ) 



(12) 
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such that the total ghost charge vanishes. Here A is worldsheet cosmological constant and 
we also have three bosonic ghost fields c{z) and c{z). We take radial ordering implicitly, 
which refiects the ordering property from noncommutative space group. 

In order to make the total superconformal ghost charge vanishing, we need vertex op- 
erators Vq in the 0-picture. We can obtain Vq by operating the following picture changing 
operator on V^i [25], 

Q = e*(e-'"<-^aZi + e2"<-^aZi), (13) 

where r]' = (1, 0, 0), r2 = (0, 1, 0) and rg = (0, 0, 1) for the components corresponding to 
the 6D compact space. Thus we have 

3 _ r ^ ^ ' 

i=i L j=i i=ijyj 

(14) 

up to the same normalization (llUj) . Here, a^^'^ is the j'th component of aj^^. Containing 
no derivatives, the higher order coupling is defined in the zero momentum limit — > 0. 
Thus the first term is not relevant. The only change is that some components of twist 
fields a are replaced by excited twist fields r, and the normalization factors (a'/2)^/^, in 
accord with the number of oscillators in flTU]) . In the next subsection, we see this change 
does not modify the calculation in the case that all the twist fields are simply a, not 
excited twist fields. 



2.2 Selection rules 

Here we briefly summarize the selection rules [3 [HI [3l [261 [27]. The vertex operator consists 
of several parts, the 4D part e*'^^, the gauge part e*^'^, the 6D twist field a'^^J, the 6D 

left-moving oscillators dZi and the bosonized fermion e**"^, as explained in the previous 
subsection. Each part has its own selection rule for allowed couplings. The selection 
rules of the 4D part and the gauge part are simple, that is, the 4D total momentum ^ k 
and the total momentum of the gauge part ^ P should be conserved. The latter rule 
is nothing but the requirement of gauge invariance. The other parts lead to non-trivial 
selection rules. In this subsection, we study the selection rule from the if-momenta and 
oscillators, as well as the selection rule from the 6D twist fields a^^J, ^) • 

The total if-momentum should be conserved like the 4D momentum and the gauge 
momentum P. For example, for 3-point couplings (V_iV_i/2Vli/2), they should satisfy 
the following condition, 

Y.r, = l. (15) 

Here we take a summation over the if-momentum for the scalar components, using the 
fact that the if-momentum of fermionic component differs by —1/2. 
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Another important symmetry is the twist symmetry of oscillators. We consider the 
following twist of oscillators, 



dZ, -> e^^'^'^dZi, dZ, ^ e^^'^'^^dZi, (16) 
dZ, e^-^^'^Z,, dZi e^2-*'^»az„ (17) 

without summation over each i = 1,2, 3. Allowed 3-point couplings (K_iVli/2Vli/2) 
should be invariant under the above Zn twist. 

However, for generic L-point couplings we have to carry out picture changing, and 
the picture changing operator Q includes non-vanishing if-momenta and right-moving 
oscillators dZj and dZi. Thus, the definition of if- momentum depends on the choice of 
the picture. However, the R-charges, which are defined as [sj^ 

R.=r,+Mi-^fi, (18) 

are invariant under picture-changing. Here we do not distinguish oscillator numbers for the 
left-movers and right-movers, because they have the same phase under twist. Indeed, 
physical states with —1 picture have vanishing oscillator number for the right-movers, 
while the oscillator number for the left-movers can be non-vanishing. Thus, hereafter A/^ 
and A/i denote the oscillator number for the left-movers, because we study the physical 
states with —1 picture from now. For simplicity, we use the notation AA/i = Mi — Mi. 
Now, the selection rule due to R-symmetry is written as 



mod Ni, (19) 



where Aj is the minimum integer satisfying Aj = 1/174, where fji = rji + m with any integer 
m. For example, for Zg-H orbifold, we have r]i = (1, 2, — 3)/6, and Ni = (6, 3, 2). 

Whereas the twist operator 0"^ itself does not transform under the twist of oscillators, 
the excited twist operator transforms like the oscillator in f|T6|) . since it is nothing but 
the product of an oscillator and a twist operator, from the transformational point of view. 
The modified if-momentum has a compensating property and the resulting amplitude is 
invariant under the twist of oscillators, as it must be because the picture changing operator 
is invariant. However the OPE dZT^ is of a similar form of a twisted operator 

9Z(^)r,(0,0)~2'=/^-V(0,0), 
aZ(^)r,(0,0)~^-^'/^i)UO,0), 

which is readily extracted from the OPE dZdZok- It has a branch structure like z^^^^'"^, 
but the part carried by from its definition to leave (l20i) . Thus the amplitude 

including excited twisted operators is the same as one including only twisted operators, 
up to the overall normalization. 



*See also [Ml and references therein. 
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We have the space group selection rule. Here, we study the selection rule for twist 
fields, 0"!^^ First of all, the product of Zn phases 7 should satisfy = 1@ Next, we 
study the space group selection rule. Now, let us consider L-point couplings of twisted 
states corresponding to {9''\Vi) {i = 1, ■ ■ ■ , L). Their couplings are allowed if the product 
of space group is the identity, i.e., 

L 

n(^'%^.)=(i,o). (21) 

i=l 

Since space group elements do not commute, nor do vertex operators, the ordering of 
vertex operators in the coupling is important. We have to take into account the fact that 
{6^ , v) is equivalent to {6^ ,v + {l — 6'^)K). Thus, the condition for allowed couplings is that 
the product of space group elements must be the identity up to such equivalence. The 
space group selection rule includes the point group selection rule, which requires Y\6^^ = 1, 
i.e. ^ fcj = (mod N) for the Zjv orbifold. The rules for the linearly combined states are 
discussed in detail in Ref [27] . 



3 The classical contribution 

Here we consider the 6D orbifolds, which can be factorized as three 2D Z^r orbifolds, 
and we concentrate ourselves to calculation of correlation functions on the 2D Z^v orb- 
ifolds. The following analysis can be extended to other cases, where 6D Z^r orbifold is 
not factorizable or 6D Zn orbifold includes 4D non-factorizable orbifold. 
The nontrivial part is the correlation function among L twist operators 

Z = (o-fciCTfca . . . Clfc^). (22) 

This can be calculated independently of the remaining components of vertex operators, 
because they commute. Some of them should be excited twisted states for the total 
ghost charge being —2. Indeed we noted above that the amplitude is the same with a 
number of factors (2a') 

From the point group selection rule we have 

E|-M, (23) 

i=l 

where M must be an integer. For the moment we assume 

M = L-2, (24) 
'"^That is automatic for physical states when the other selection rules are satisfied ^51 [?7] . 
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and we will relax this condition later. This choice is the most convenient one because 
of two reasons. First, by doubling trick we can relate the corresponding amplitude with 
that of open strings [T7j, where this is the closedness condition for L sided polygon. This 
gives rise to the generalized Schwarz-Christoffel transformation and the area rule that we 
will show below. Another reason is that we can obtain the closed form of integration, in 
terms of a multivariable hypergeometric function [29J. 

If we use path integral formulation, the correlation function fl22|) is divided as 



2 = 2qu ■ ^ exp(-S'ci), (25) 

{Zcl} 

according the classical value and the quantum fluctuation around it, i.e. Z = Zc\ + Zqu- 
The corresponding tree-level Feynman diagram is a sphere with a number of vertex 
operators inserted. Every object in string theory, including the vertex operators and the 
correlation function (!22l) can be separated into holomorphic and antiholomorphic part. 
Considering one of them, holomorphicity restricts many things in a very simple form. 
Using the compactification C U {cxo} = 5*^ we can cover all the coordinate on the sphere 
by holomorphic coordinate z except infinity. To take care of infinity we introduce another 
patch z = 1/u where z ^ oo becomes u = 0. Considering a holomorphic solution from 

d^Z = -z^dZ (26) 

if we need the LHS well-behaved at u = 0, in the RHS dZ should drop faster than z""^ as 
z goes to infinity. 

3.1 L-point coupling 

We begin calculating the classical contribution first. We see that the classical solution 
is the completely factorized part for each inserted operator. Therefore, from the OPEs 
([3]) and (jl]), the classical solutions are obtained by the holomorphicity and the desired 
singular structures as 

dZci{z) = auj{z), 

L-2 ^27) 
dZa{^ = Y,h^'\^)- 

1=2 

Here we define the basis of [L — 2) functions 

L 



\ki/N-l 



UJ[Z] = ^^[Z - Z, 
i=l 

L L-2 



(28) 

n(^-^.)-'^/^ n (^-%)' ^=2,...,L-2. 



i=l j='i,j¥=i 
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For the symmetric orbifold they are complete since dZci{z) = {dZci{z))* and dZci{z) = 
{dZc\{z))*. From (l26l) . the whole part fl27|) should behave as z'"^ as z oo. Using 
(!24|) . we can see u^z) does. For lj"-{z), we admitted an additional degree of freedom, i.e. 
changing the power of z — zj singularity by integer, since it does not modify the branch 
cut structure. Hence, we have many free parameters 6;. 

To determine the coefficients a, bi, we should consider the global monodromy condition. 
The relation ([2]) does not take into account the global phase if we transport branch cuts 
from more than one fields. Taking a contour C encircling more than one points gives the 
relation between overall coefficients in ( !27l) and net translation v in the target space 

AcZ = i dzdZ + i dzBZ = v. (29) 
Jc Jc 

This expression makes sense only if there is no additional phase. The quantum part does 
not carry any amplitude AcZq^^ = 0. Hence, Eq. ( l29l) shows a purely classical contribution. 
Upon integration relating these branch cuts, the Pochhammer loop [30] is a clever way 
to encompass the two branch points nontrivially without phase. For each branch cut the 
contour goes in and out exactly once through the cuts, depicted in Fig. [H Its effect is 
to encircle the fixed points: fi clockwise, /2 counterclockwise, fi counterclockwise and 
then /2 clockwise. In terms of space group elements {uj,v[ = (1 — u;)(/i + vi)) and 
{9, v'2 = {I — 0){f2 + V2)) with f 1 2 G A, we obtain 

ie,v'2)iiu,v[)-\e,v'2r\uj,v[) 

= {l,{l-uj-')v2 + {l-0)oo-hi) (30) 
= {l,{l-u;-'Kl-9){f2-h + v)), 

where v = V2 — Vi. The encircling is not necessarily once, i.e, u ^ 6 in general, in which 
we cannot draw branch cuts. The net effect is pure translation. It turns out that every 
contour is generated by the basis of Pochhammer loops Cj encircling i-th and {i + l)-th 
points. 

Taking Pochhammer loops Ci around the vertices Zi and Zj+i, from (I5UI) we obtain 
Ac.^ci = / dzdZdiz) + (f dzdZdiz) 

^ _ g-2..fc,/JV)(^ _ g2..fc,+,/JV)(^^^^ -f^ + v) (31) 

= 4e-('^-'=-^)/^ sin sin (/.^^ -f, + v). 

We have (L — 2) vectors /j+i — fi + v and L angles (with the constraint fl2^ ). which 
completely specify L-sided polygon. 

Later we can express the solution in terms of the following integrals 

= (f dzuj{z), Wi = (f dztuiiz), (32) 
Jci Jci 
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Figure 1: Pochhammer loop. Whatever direction the branch cut we have, the nontrivial 
loop goes in and out respectively exactly once through each cut, enclosing two branch 
points. 



and 

Wi = (1 - e-2-*/^)(l - e2-^*^«+i/^)F/, i = 1, . . . , L - 2, (33) 

with 

Fl = / u{z)dz, Fl = / Lj\z)dz, i = 2, . . . , L - 2. (34) 

Note that F^ and form (L — 2) x (L — 2) matrices. In Appendix, they are expressed 
in terms of multi- valued hypergeometric functions [T71 121]. With SL{2, C), we can set 
^ly zl-i, zl to be 0, 1, cxo respectively and the others to the cross-ratios of Xj. 
Plugging (l27Il . the solution is expressed as 

CiFi = fi+,-fi + v, Z = l,...,L-2, (35) 

where we defined Ci = a,Ci = b*, and by inverting them we obtain 

L-2 

c, = Y,Ui+i-fi + v){F-% (36) 
1=1 

where the inverse is taken with respect to the matrix basis with indices /, i. Plugging into 
the classical action, we obtain the final solution 



sM„...,x,.,) = ^(^a\^i+j2b:b,n^ 



(37) 



where 

/(X2, . . . ,Xl-2) = / (i^z\uj{z) 



I',j{x2, . . .,XL-2) = / 

Jc 



c 

'2 -i 

a 

c 



(3^ 

zuj\z){uj' {z)y . 



10 



We can expand this action by products of holomorphic and antiholomorphic functions, 
with careful choices of contours. It is nothing but the relation between open and closed 
string amplitudes before integration over XiS [3T] . 



where 



L-1 

1=2 
L-2 j-1 

j=2 i=0 



L-1 



k, 



1 - exp(-2m 2^ — ^ 



1=2 



1 — exp(— 27ri 



i=i+l 



(39) 



p.p. 



-oo .-I 



Fl-1 



2 1 
OO 



L-1 



dz\{[ 



z — z 



1 

-2, 

-{l-k,/N) 



(40) 



Plugging in (j37l) . we obtain the classical action 5'ci(x2, . . . ,xl-2)- It it a function of 
L — 3 complex variables Xi which will be integrated out in the final amplitude. Later, we 
will integrate this with variables X2i ■ ■ ■Xl-2 over the entire complex plane. Among these, 
using the saddle point approximation by adjusting XjS or equivalently ratios = 
2, . . . , L — 2, we find a minimum 



Fj+l _ fi+2 — fi+l 

Fi fi+i — fi 



(41) 



Note that in general the integrals are complex and we coordinated the fixed points as 
complex vectors on a given 2D torus and orbifold. 

Inserting these into (136!) . we have q = 6* = for all i > 1, except a = ci 7^ 0. The 
solution is nothing but the generalized Schwarz-Christoffel transformation [29], whose 
original version maps the upper half plane into inside an L-polygon 



dZc\{z) = aLj{z). 



(42) 



Namely, the points Xi are mapped to vertex Zi and the turning around angle is given by 
kiir/N. In this case, we obtain the instanton contribution is exponential of the polygon 
area 

-("area of the polygon" 



exp 



2'Ka' 



(43) 



This is valid under the assumption (^^, i.e. forming a polygon, but in general case we 
have more fundamental interpretation shortly. 
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Of course, there are other minima with the same value, where a and all bi vanish 
except one, say bk ^ 0. This corresponds to F^^^/ Fl = (/j+a - - fi) and the 

Schwarz-Christoffel transformation corresponds to dZc\{z) = blu'{z)^. 

In forming the area from the Schwarz-Christoffel transformation, the ordering is im- 
portant. If we just exchange two fields, we cannot satisfy the space group selection rule in 
general, and the polygon becomes self-crossing, where the area rule is not applicable. In 
the correlation function we take the radial ordering. In the superpotential of effective field 
theory, we do not see the ordering, since the integration over all Zi completely symmetrize 
the amplitude. 



3.2 Four-point correlation function 

The four-point correlation function provides a good example of calculation of the classical 
part. In this case, the functions Fi are well-known hypergeometric functions, which are 
solutions of second order linear differential equation. It is known [32] that if any three 
of the solutions have the common domain of existence, there be a linear relation among 
them. In our case, we can express all of Fi in terms of, say, Fi and F2. They are shown 
in (UnHD and ffT09D of Appendix. 

Plugging these to ( l39ll the holomorphic part is obtained as 



where 



Cll 



I{x) = cii|Fi|2 + C12F1F2 + CI2F1F2 + C22IF2I2, (44) 

sin(7rfci/A^) sin(7r(fc2 + fcg) /N) 
sin(7rfc4/iV) ' 

sin(7r/c3/iV) sin(7r(A;i + k2)/N) 

C22 = ■ . , , 45 

sm(7rA;4/A') 



C12 = -e^*^^/^ 



■. , , , sm{n{ki + k2)/N)sm{n{k2 + h)/N) 
sm[7rk2/i \ ) H 



sin(7rfc4/A^) 
The coefficient C12 reduces to 

= e"^^/^ sm{TTki/N) sinink^/N), (46) 

only for the polygon case, using (IMl) . The prefactor e^'^^'^l^ in C12 is the relative phase of 
(complex numbers) W32 and V21, where % = fi — fj + v with t> G A. From (!36i) we have 
the coefficients 

V32F1 + V21F2 



-I 5 



F1F2 — F2F1 ^^Y) 

^^32-^1 - ^^21-^2 

i^li^2 ~ -^2-^1 



12 



We can obtain the antiholomorphic action /'(x) and integral Fl{x) from I{x) and Fi{x) 
respectively, by substituting ki/N 1 — k-jN and x ^ x. With these we obtain the 
classical action (1371) . The action does not have manifest duality symmetry, since we have 
fixed four points by SL{2, C). 

We define the following modulus 



r 



El 



(4J 



which is in the Z2 case the modular parameter of two-torus, made by connecting two 
Riemann sheets with two branch cuts [6l[7j. As expected from fHTI) . the minimum of I{x) 
is obtained at 

VZ2 



T 



V2I 



Thus we have a 



S, 



-V21, b = and the minimum of classical action is obtained as 
1 



cl,min 



2na' 



Cll|t^2lP + £12^21^32 + CI2V21V32 + C22|^^32P 



(49) 



(50) 



For the case of polygon, the classical action Sd^min reduces to 



27ra' 
1 

"27ra' 



|t;i4p sin(7rfci/A^) sin(7rA;4/iV) |t;32p sin(7rA;2/A^) sin(7rA;3/A^) 

2 sin(7r(A;2 + k3)/N) 
\v2i\'^ sm{TTki/N) sm{iTk2/N) 



2 sin(7r(A;i + A;4)/A^) 
|i;43p sin(7rfc3/A^) sin(7rA;4/A^) 



sin(7r(A;3 + ki)/N) 



sin(7r(A;i + k2)/N) 



(51) 



This is the area of the quadrilateral formed by vertices at the fixed points /i, /2, /s and 

u 

In the case with ki + k^ = N and/or ^2 + ^3 = N, this expression is not well- 
defined. Without loss of generality, the case with ki + ki = N and k2 + k^ = N leads to 
ki = ks = N — k, k2 = k4 = k, by use of flMl) . and such a case has been calculated in [?]. 
Here we have to come back to (jSj), and the result agrees]^ The case with ki + k^ = N 
or k2 + k^ = N leads to ki = N — k, k2 = k, k^ = N — I, k^ = I, and such a case has been 
calculated in [H]. 

From four-point amplitude, we can obtain three point amplitude by taking x to, say, 
00. In this case the fixed points /3 and become coalescent and the classical action 
'S'ci,min reduces to 

1 |^32psin(7rA;2/A^)sin(7rfc3/Ar) 
"^cLmin 2^^, 2 sin(7r(A;2 + A;3)/iV) ' ^ ' 

Note that this action depends on the choice of contour "picture". Here, we chose one 
encircling two fixed points and f2- We do not need worry about whether v is actually 



To compare between our results and 7], we have to replace our modulus r by e'^*'-'^/^ ^V. 
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compatible to factorization The Pochhammer loops in which v^2 and fi4 belong are 
independent. 

In the special case with ki = k2 = = = N/2, we have cn = C22 = 0, C12 = —i and 
Fi' = Fi,F^ = -F2 yielding 

-2t;2i— = - — -\v2iv-i2\, (53) 



27ra' z 27rQ;' 

which is again interpreted as twice the area of the rectangle, in unit of 27ra' if V21 is 
orthogonal to v^2- This is the case of order 2 subsector (A^/2-th twisted sector) in even 
order orbifold. Note that this action is not the minimum action, since in this case the 
classical action is not the function of r. In this case the area rule interpretation is somehow 
ambiguous. We will study more detail in the following subsection. 

We have considered L-point couplings only for L = 2,3,4 as examples. However, 
we will study that higher order L-couplings reduce to a combination of lower order L'- 
couplings with L' < L hj the discussion of field coalesce in subsection 3.4. In addition, 
since the number of fixed points on T'^/Zn orbifolds is limited, we can expect that most 
of higher order L-couplings can be written as combinations of L-point couplings only with 
L = 2,3,4. We will study this expectation in separated papers ^27^ |33j, by examining 
concrete orbifold models. 

3.3 Non-polygon case: the meaning of area 

We assumed the polygon condition (p4l) is satisfied. However, in general, the following 
relation 

E|si-2, (54) 

1=1 

is possible. In the inequality case, the holomorphic part of the classical solution decays 
faster than z~'^, whereas the antiholomorphic part decays not faster than The only 
sensible way of treating is to make antiholomorphic part vanishing. 

For example, in the Z4 orbifold, the coupling of four first twisted sector fields, 

shown in Fig. [2]^a), satisfies the space group selection rule. Note that o'{efi)(^{0,o)(^(e,ei)O'{e,ei) 
does not satisfy the space group selection rule, and the ordering is important. All of twist 
fields are of order four, which cannot satisfy the relation fl2^ . The classical solution is 
obtained as 

aZei = a{z - z,r'^\z - Z2)-'/\z - z^r^l\z - z,)~"\ 
dZ,i = 0. 
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2,4 




1,3 



(a) 





(c) 



Figure 2: (a) On Z4 orbifold, four order 4 couplings: two at a vertex, and the other two at 
the opposite vertex, (b) Massless strings are locahzed (c) Quantum effect grows twisted 
strings to form a untwisted string, which is a intermediate state described by instanton. 
The shaded area, where V21 is the diagonal, is the minimal area swept. 







Figure 3: Naive area rule is the special case of the swept-are rule when the twisted strings 
form a polygon. 



From global monodromy condition, we have 

and the classical action is given with cn = C22 = 1, C12 = — 1(1 + i), yielding 

1 = 2. 

Thus we have the classical contribution 



5, 



1 



cl 



r2k2i| 



27ra' ' 

which is interpreted as the area of the square whose diagonal is V2i- 
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Similarly, for the coupling 



O"(e,0)O-{e2,ei)O"(e,ei), 



we obtain the classical action, 



In these cases, we have a different area rule: The "area" is not that surrounded by 
fixed points, but one as follows. The classical solution describes a local minimum of 
the action, which is the instanton of worldsheet nature, suppressed by a'. The selection 
rule tells us that these twisted strings can potentially make an untwisted string, which is 
not possible due to energetics for massless strings, since they are completely localized at 
certain fixed points, as in Fig. [21(b). However they can oscillate to grow to be large size, 
and above a certain threshold, they can form an untwisted string as in Fig[2]^c). Noting 
that the instanton describes tunneling between vacua which is energetically forbidden, we 
can understand that forming untwisted string corresponds to such tunneling. 

Still we can have the hint from the modified area. It is the sweeping area for localized 
twisted strings to grow to become a untwisted string. For the polygon case, i.e. that 
satisfying the condition fl2^ . this interpretation is still valid, since still the area swept by 
twisted strings at each vertex makes the polygon area. 

On the other hand, in the case not satisfying the condition fl2^ . we lost the interpre- 
tation of the mapping dZ being a generalized Schwarz-Christoffel transformation, since 
in the target space the fixed points fail to make a polygon. 

The subleading correction is generated by identical fixed points on the orbifold, but 
more separated in the covering space, i.e. / + (1 — 6^)K. Since they are identical points, 
they satisfy the selection rule. With this interpretation, we can understand the classi- 
cal solution which makes more than one untwisted strings possible. For example, for 
<^(6i,o)<^(e2,ei)<7(e,ei), the coupling cr(0 o)C"(e2,3ei)<7(e,3ei)5 where each fixed point belongs to the 
same conjugacy class as the previous one, satisfies the space group selection rule, but 
corresponds to a large V21 and a large instanton action. 

3.4 Fields coalescent at the same fixed point 

In most cases, some of the fields sit at the same fixed point /. The corresponding corre- 
lation function might be obtained by taking the limit Zj Zi in the correlation function. 



Note that this limit is not always well-defined. It is because twist fields do not commute. 
By conformal symmetry, the OPE has the generic form. 




(56) 




(57) 
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Equating the conformal weights of the both sides, we have 



— h — h — /) —J TVTV — (<^R') 

'^u- - - _!)(!_ I) (A;, + A;,- > iV) ■ ^ ^ 

Because of the nontrivial branch cut, this relation is asymmetric under the exchange of 
two twist fields. This property is also reflected in the space group elements, which do not 
commute, either. We can define an invariant block of twist fields, which correspond to 
the identity of the space group (1,0). Thus, these invariant blocks commute and satisfy 



E 



Kij = mod 1. (59) 



In the case where two points Zi and Zi^i are in the successive order, and we can merge 
two twists without ruining the radial ordering. They are neighboring points as polygon 
vertices. Then, from fl57p the two twists reduce to a single twist with the summed order. 
Also this implies the classical solution becomes 

■■■{Z- Zi)''^/^-\z - Zi+i)^^+'/^-^ >---{z- ^^)(fc»+fc«+i)/^-2 . . . (gQ) 

where, again, the even integral power is not relevant to branch structure, so that we can 
make ki/N + ki+i/N — 1 lie in [—1, 1) for instance. This means, not all of L vertices form 
the L-polygon, but effectively one with the lesser vertices. (Recall that the positions of 
vertices are given as singularities in the classical solution.) In fact this is the familiar 
case when we obtain a three point function from the four point function by setting two of 
the points coalesce. In the latter limit, the polygon is triangle. We can see this in terms 
of space group elements. Neglecting gauge group, which is not involved in the classical 
amplitude, we cannot distinguish the product 

(to , vY and {uj'^ , V + ujv) , (61) 

when two identical twisted fields sit at the same point, and are put on the neighboring 
points in the correlation function. 

There is the case where ki + fcj+i is integer. In this case there is no singularity in the 
classical solution and such vertex does not contribute to the area rule. In the extreme 
case where k order-/c twisted fields sit at the same points, the coefficients of higher order 
couplings are not suppressed. Thus all we need to consider is the other nontrivial cou- 
plings. Fortunately, not all of them survive: From i?-symmetry invariance, the correlation 
function is further constrained. 

However, from the radial ordering, there is a case in which the exchange of two branch 
cut fails to give well-defined radial ordering. This is not possible for two fields which 
are not the neighboring fields in the correlation function, since the space group and the 
mapping of classical solution is not commutative. 
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For example, in couplings among twist fields corresponding to two {6, 0)'s, two {6, ei)'s 
and two {9,ei + 62) 's on the Z3 orbifold, there are two following combinations possible, 
satisfying the space group selection rule 

0" (e,0)O"(e,0)O" (6»,ei)0" (6l,ei)Cr(6l,ei+e2)0"{6»,ei+e2) , (62) 
{e,0)O"(e,ei)O"(6l,ei+e2)'^(6»,0)Cr(6l,ei)Cr(e,ei+e2) • (63) 

In the first case, the solution behaves like 

dZ = a{z - zi)-^'\z - Z2)-^'\z - z-i)'^l\z - z^)-^'\z - z,,)-^'\z - z^Y^'K (64) 

The two points, zi and Zj+i for i = 1, 3, 5, correspond to the same fixed points on the 
target space, i.e. Z{zi) = Z{z2) = /i, Z^z^) = Z^z^) = /a and ^(zg) = Z^zq) = f^. Thus, 
we take the limit Zi^i zi for i = 1, 3, 5. In such a limit, the solution behaves like 

dZ = a{z - z^)~^'\z - z:i)-^l\z - z^)-^'\ (65) 

Thus, in the space group point of view, it is indistinguishable from the coupling among 
the second twisted couplings 

O"(e2,0)O"(e2,ei+e2)O"(e2,e2), (66) 

which gives only the area of triangle. On the other hand, the latter coupling fl63p gives 
twice the area of the triangle. The classical solution maps from the different points to the 
same points. For example the holomorphic part behaves 

dZ = a{z - z,)~^'\z - z^r^l^z - z^Y''l\z - z.r^l^z - z,)-^'\z - z,)-^l\ (67) 

and the two points, Zi and Zi+^ for i = 1,2,3, correspond to the same fixed points on 
the target space, i.e. Z{zi) = Z{zi) = fi, Z{z2) = Z{z^) = /a and ^(zg) = Z{zq) = /g. 



18 



Otherwise the ordering is ruined. It is also understood that two different and independent 
sets of twisted fields can sweep the triangle, which means the instanton corrections from 
the complete polygons are additive in the action. This case is depicted in Fig. HI 

We summarize how to calculate the classical part of L-point couplings among twist 
fields (Tfc. {i = 1, ■ ■ ■ , L). First, we classify all of possible ordering of these twist fields, 
which satisfy the space group selection rule. For each ordering of twist fields allowed by 
the space group selection rule, we consider the following procedure. We combine two or 
more twist fields sitting at the same fixed point to a single twist field like fl57|) and Fig. 3 
if possible, that is, they satisfy ( !59|) . When their total twist is just (27m) twist with integer 
n, correlation function reduces to much simpler form. The resultant correlation function 
can be written as product of invariant blocks, which satisfy the space group selection 
rule like Fig. 3. Each block includes smaller number L' of twist fields|l| Then, for each 
block, we calculate classical contributions, e~'^'=', i.e. instanton actions corresponding to 
the minimum action and larger ones, and take their summation, i.e. J2i {vi+i j 6xp(— Sd). 
Next, we take a production of classical contributions corresponding to each block, i.e. 

n I E exp(-5e,)J. (68) 

invariant blocks \ , } / 

Finally, we sum over all of possible ordering to obtain the total coupling, i.e. 



E 

possible ordering 



n E ^M-Sci) 

invariant blocks \ i.jiii^i ^} 



(69) 



Because of the ordering, the correlation function does not possess the worldsheet 
duality like in the Virasoro-Shapiro amplitude [31] . In the vanishing momentum limit we 
do not distinguish the channel, thus the effective coupling like Yukawa coupling does not 
distinguish the order. However even in this case, the above two cases are distinguishable, 
since this contribution is the worldsheet effect, suppressed by a'. 



3.5 Linearly combined states 

In higher order twisted sector of a non-prime order orbifold, there are states formed by 
linear combinations, as in ([6j) , due to the orbifold projection. The linear combination 
of states can be more precisely defined by that of vertex operators. It follows that the 
corresponding classical solution consists of linear combination of individual solutions be- 
fore combination, and also the relative weights are inherited. For example the classical 
solution involving ([6]) contains the factor 

^ ((. - zi)-'^'/'' + l{z - .2)-^+'=/^ + • • ■ + l'-\z - z,r'^'"^^ , (70) 
^ In most of cases, L' may be equal to L' = 2, 3, 4. We will study this point in concrete models [55] . 
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where Zi is mapped to the fixed points {9, 9^ ^v). 
For such couphngs, we observe two points 

1. The selection rule for a linearly combined state is derived from that for each term. 

2. A part of the classical solution, which does not satisfy the selection rule vanish. 

The proof for the first is given in Ref. p7j, and we easily see the latter is the case. We can 
show that such term not satisfying the classical solution always contain the same element 
more than once. Setting these to the other(s) corresponding to shrinking the area to zero, 
if we do not want to change the angles corresponding other vertices. 

For instance, in the second twisted sector of Z4 orbifold, we have a coupling including 
linearly combined states 

^(e^fl) ■ (^(e2,ei) ± 0-(e2,e2)) " ^(02,ei+e2) " (^(e2,ei) ± {9^ ,e2)) ■ (71) 

The total coupling is given by summation over each of the four correlation functions. 
However, the nonvanishing contributions come from the only ones satisfying the selec- 
tion rule. One that does not satisfy the rule, e.g. (T(6i2,o)C'"(e2,ei)O"(02,ei+e2)<7(e2,ei), does not 
contribute. 

To sum up, for the linearly combined state, we can treat a linearly combined state 
as a complete physical field, the classical amplitude has only contribution from the parts 
satisfying the space group selection rule. 



4 The quantum amplitude 

Now we determine the quantum part of the amplitude fl25|) . 



(72) 



We use the stress tensor method [3 [HI [IZ], in which (1721) can be indirectly calculated, 
relying on only the holomorphicity. All the information can be read from the Green's 
function 



{-^d^Zd^ZYlak^) 



qu 



Its form we know from the holomorphic part of the OPE ([3]), 

{z-Xi){z- Xj) Uk^iji^ - Xk) 



g{z,w) = u!{z)uj'{w) 



.i<j 



[z — w) 



(73) 



(74) 



where uj'(w) = Yli^ ~ ^i) ^^^^ ■ The condition i < j avoids overcounting, although 
otherwise the formula would be more symmetric. The prefactor uj{z)uj'{w), from (l28l) . 
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gives the desired pole structure. The normahzation is determined by requiring the z w 
hmit of g{z,w) to be 1, and the condition 



(75) 



leads the desired conformal weight hk, as the coefficient of each double pole in z — w. To 
have no residue in z — Zi, we require the condition for aij, 



L 



1 — - 

iV' 



(76) 



where we use the fake number a^j = aji for the case i > j. Summing over i, this condition 
implies (175|) . We have less condition for many aij, reflecting the freedom of our choice of 
A. Of course, the physics does not depend on the choice aij. As z ^ w, we have the OPE 



qu 



IE 



kj 



2j^NNiz-z,)iz-z,) 



2 ^ 

Kj 



1 1 
+ 



2' Zj 2' Z rj 



+ 



A 



(77) 



Uiiz - Zi 



We extract the residue {z — Zk) for fixed k, and take z ^ Zk limit. 

In the limit w ^ z, the Green's function becomes energy-momentum (stress) tensor 



Tiz) 



~ lim 

w—*z 



-\d.zd^z- ^ 

2 [z — wY 



(7J 



where the last term arises from the normal ordering. Sandwiching the OPE in the corre- 
lation function 

T(z)afc,(zfc) ~ — + , (79) 

{Z - ZkY Z- Zk 

with the fixed ak^Zk) and its conformal weight hk given in Eq. ([5]), we can completely 
calculate the holomorphic part up to normalization 



d 

dz. In = = lim {z - Zk) 



We obtain 



i<j 



-qu 



an — 1 — 



{nz)u^ k.) 

crh 

^qu 



qu 



hi 



ki 
N 



I — hi 

N 



+ 



{z - ZkY 



A 



^0) 



Zk - Zj ]\i^k(Zk - z, 



(81) 
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At this stage, using SL{2,C) symmetry, we set zi ^ 0,Z2 ^ x,zl-i ^ 1,zl ^ oo. 
Then zl dependent terms vanish since zl — ^ oo, but there are corrections to this from 
A, which are also dependent on zl- We now have L unknown for L — 1 constraint, and 
remaining freedom corresponds to redefinition of A. 

Similarly we define the antiholomorphic Green's function 

L-2 

h{z,w) = J2 Bijuj''{z)uj'.{w). (82) 
To obtain A and Bij, we apply the global monodromy condition for the quantum part 



c, 



dzg{z,w) + (p dzh{z,w) = 0. 



(83) 



Ci 



It is convenient to call the first term in the RHS ([71]) as gs{z, w). If all the Zi are different, 
the matrix has the inverse. Thus we can multiply it to eliminate A and Bij, so that 
we obtain 

L~2 



u;'{w)A = -J2iW~'yi / dygsiy 



L-2 



L-2 



Y,ooHw)B,, = -J2{W-')\ I dygs{y,w). 
j=i I -^(^i 

Since we set z^ — > oo, the only relevant terms are ones containing the factor aiL{z — Zi). 
We divide them by uo'{w) and take the limit w — >■ oo. Then, in the sense of (!79l) . we read 
off the residue of {z — z^) in the integrand 



hm {z - Zk) 



gs{z,w) 



lim 

w^oo Uj'[W) 



a,- 



ki 



N I Zk - Za 



^5) 



where we have replaced with by the relation (1761) . Therefore we have no dependence 
on the specific choice of gs{z,w), or aij, as it should be. 
Integration around the contour Ci gives 



+ dz. In 



Ylizk - Zj 



Zi Z'j 



(86) 



The first term in the RHS is contained in the derivative of det W = det Wl'. By the chain 
rule, we obtain 



L L-3 L 

j=l,j^k 1=1 



^7) 



1=1 
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For the second term in the RHS of (!86|) . in the same way, we find the singular structure 
as z ^ Zk 5 

Thus we obtain 

L-3 L 

n ~ n ~ 



a,,lnZj, = 9,,(detl^)-Mn 



.ij=l,i<j j<j=l 



9) 



up to antiholomorphic action. In fact, there are additional factors, since we have used 
only the part of the terms Ylii^'^Yi^Zk^l ( |871) . The remaining terms completely 
vanish if we calculate the antiholomorphic part in the same way. The similar analysis 
is carried out for Bij, by dividing by uj''^{w). Combining it with the holomorphic part 
^qu = -S^Ju ■ ^qn'^^ we arrive 



L-2 L 



Z^,= {detW)-' n (^*-%) n {zi-z,)-^'^^^^'-^\z.--z,)-^^, (90) 

i,j=l,i<j i,j=l,i<j 

up to the normalization, to be determined in the following section. The various branch 
cuts refiect the noncommutativity of the vertex operators under permutation, due to 
orbifold phase. As well as the classical part, the quantum part has no world-sheet duality, 
either. Nonetheless the complete amplitude will be single-valued, completed with other 
parts of vertex operators. 

5 Factorization and normalization 

To normalize the amplitude, we choose the reference of normalization 

{akiz,z)aN_kiO,0)) = 1 ■ \z\-^^^'-^\ z ^ 0. (91) 

Equivalently, two twist fields of the opposite twists coalescent at the same point become 
the identity operator. In general we cannot make use of the normalization fl9Tl) . since for 
each twist field we need one with the opposite twist. Thus we use a doubling trick. 
First we calculate 2(L — 1) point function with twists 

5(1; 2; . . . ; 2L - 2) = {ak^ {zi)aN-ki (wi) . . . crk^_^ {zL-i)crN-kL_, (^L-i))- (92) 

By setting wl-i — ^ zl-i we obtain the product of a (2L — 3)-point function and a 3-point 
function. From the unitarity, the intermediate state gives an on-shell pole, with the lowest 
state being massless gauge bosons, 

g(l;2;...;2L-2)^./^ ^ft^^-^^^-^;:^)f^^^-^-^^-^' . (93) 

J [271)^ —k'^ + ie 
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Figure 5: Factorization of 2(L — l)-point amplitude to [2(L — 1) — A;]-point and fc-point 
amplitudes. Unitarity tells us the intermediate states have mass poles, whose residues are 
the product of the factor amplitudes. The correct normalization does not depend on a 
specific factorization and choice of intermediate state. 



This factorization is schematically drawn in Fig. [51 

Concentrating only twist fields, the intermediate gauge boson does not contain a twist 
field. The former contains only 2L — 4 twist operators and the latter has the form as in 
( !9T|) . Therefore we have 

{ak^aN-k^ ■ ■ ■ ak^_^aN-kL-i) = {cTk^cFN-k^ ■ ■ ■ aki^_^aN-kL-2) {(^H-if^N-kL^i) ■ (94) 

We have already known the general solution. For this we have the same action, including 
the classical and the quantum parts, except the factor det W . In this, each term Wl^ 
contains the common factor 

1=1 ^ ^ 
Thus there is a discrepancy between the two by the factor 

= 4,.(..-..-0V- ^ sin (96) 
detiy(/_i) N N ' ^ ^ 

where the subscript in the determinant indicates that these are for / and / — 1 point 
correlation functions, respectively. This is so, since the normalization of 2 (L — 1) point 
function have no special dependence of a specific fcj. In effect, we replace det W in 
fl90|) with det F, which is defined in (^^. This is a desirable result, since any choice of 
factorization should give the same result, not depending on the special set of contour 
choice, which is reflected in sinoj. 

We cannot go from the four-point function into the product of two two-point functions, 
since the point z cannot be arbitrary close to all of 0, 1, oo, already fixed by SL{2, C), at 
the same time. For this, we do the Poisson resummation |7j. Since it involves a lattice 
transformation into its dual lattice, as a result we have the overall factor of the lattice 
volume. 
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Finally we set zl, -Zl+i, • • • ^2L-2 to infinity to obtain the OPE 

(6»'=,t;)eP,sel 

where the selection rules should be satisfied. The RHS is nothing but the product of two 
identical (L — l)th order couplings. Comparing the coefficients, we obtain the coupling 
of order L 

^L-point(/si, /^2, • • • , ki) = lim [det F^^2L-2)ih, k2, . . . , ki-i, N-ki, N-k2, . . . , N-kL-i)] 

all Wi^oo 

(98) 

where 

^=L-l-y^. (99) 
i=i 

For the rest of vertex operator components, we have 

{c{Zi)c{Zi)c{z2)c{z2)c{z3)c{z3)) = \Zi - Z2\^\Zi - Z^\'^\z2 - 2:3^, 

(e-^zi)e-^^^2)e-^^Z3)) = (% - Z2)-"\zi - Zz)-"\z2 - ZzY^". 

{riz2Wizs)) = r''iz2-zsr\ 

L 

'J[e''-''{zi,z,))= n \zi-zjr''-'^, (100) 



i,j=^,i<j 
L 



ii'-^-''iz.))= n ( 

«j=i,«<i 



z, - z.T''^-'^/^. 



Multiplying these, and using the massless condition, the branch cuts disappear in the 
overall amplitude. Including universal geometric factor of sphere 8'n'g~'^a'~^, we have the 
overall factor 

2(47r)-^+Vo-2e(^-2)*a'l(^-4), (101) 

up to contributions from picture changing and/or oscillator excitations (fTOj) . For the four 
point coupling, we have volume factor suppression from the compactification. 

As the simple example, we can extract the information for three point correlation 
function. By SL{2,C) symmetry, we can fix three arbitrary given points completely, 
thus the correlation function contains no information. To obtain it, we need four point 
correlation function with twists 

(crfc(0, 0), aN-k{x, x), a/(l, 1), o-Ar_,(oo, oo)), (102) 
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where we consider only one two-torus. From (1321) . we have 



det F{x, x) =FlF^ - F^F^ 

=i?(l-|,^)F(] 

— {k ^ I, z ^ z exchanged) . 



k I . -I I I 



-• 1 - z)F(— 1 



Af ' -^i ^/ 



(103) 



Let us assume k + I < N. In the hmit z, z ^ oo, using the relations flllOp and fillip in 
Appendix, the most dominant part has the coefficient 



det F ^ TT 



r2(i-|)r2(i-±) 



-ik-K 



n -1 



sm 



Td 

N 



+ 



sm 



TV 



(104) 



where we used the relation r(a)r(l — a) = vr/ sinavr. Thus, taking the entire compact 
dimension, the Yukawa coupling is obtained as 



d/2 



k,l 

fajbjc 



5'YM,4D 



n(2A 



Aj 



,1/2; 



r(i-|)r(i- 



r(i 



N 



sm 



■Klj_ 

N 



+ 



1 1/2 



sm 



J 



E 



exp( 



-5, 



cl 



(105) 

with the classical action given in fl52l) . For the case fcj + Ij > N, we can obtain the 
corresponding amplitude by replacing kj, Ij hj N — kj, N — Ij. It is notable that, in the 
quantum amplitude, there is no contribution from the geometry, since it is canceled by 
the same dependence in the four dimensional gauge coupling arising from the dimensional 
reduction. 

We can obtain similar expressions for the higher order couplings, expressed in terms 
of multivalued hypergeometric functions. One can be convinced that the asymptotic form 
of generalized hypergeometric function in the above limit is the ratios of gamma functions 
[35] with arguments of C(O.l). Thus we expect that the quantum parts of the higher order 
couplings are roughly of order one. Thus we see that the size of higher order coupling is 
dominated by the classical part. 



6 Conclusions 

We have calculated couplings of arbitrary order, among untwisted and twisted fields of 
heterotic string on orbifolds, using conformal field theory. In the low energy limit, they 
correspond to the matter superpotential. They are given by the zero external momentum 
limit of radially ordered correlation functions. The specification of orbifold and the shift 
vector determines the possible couplings. This provides us with lessons for constructing 
low-energy effective field theory, in particular for vacuum configurations of compactified 
string theory and realistic quark and lepton masses. 

The higher order couplings are complicated due to two things. The ffist is the technical 
difficulties dealing with arbitrary number of twist fields and even some excited twist fields 
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from the picture changing. In the calculation of higher order coupling, the latter has 
just an effect of changing the normalization, not changing the transformation property 
and the branch cut structure. The other difficulty arises from linearly combined states, 
which appear in higher twisted sectors of non-prime order orbifolds. They make the 
interpretation of the classical and the quantum somewhat tricky. 

The selection rules (mainly studied in Ref. [27]) from the location of fields and the 
i?-symmetry can be possible origins of discrete quantum numbers in effective field theory. 
They also provides the understanding on discrete flavor symmetries. Because of the 
restrictive form of selection rules, in particular in the heterotic string models, only limited 
number of couplings are possible, since there are limited number of orbifolds and fixed 
points. This will be classified elsewhere [33] . 

The classical part is an instanton amplitude of world-sheet nature. Its size is expo- 
nentially suppressed by the effective area swept by the twisted strings to form untwisted 
strings. The latter is energetically not allowed, metastable intermediate state. This gener- 
alizes the naive area rule even if the twisted strings do not form a polygon. Decomposing 
the locations of twisted fields, we can obtain handy rule for calculating the size. For 
couplings involving linearly combined states, the only contribution comes from the terms 
satisfying the space group selection rule. The others vanish individually. 

We have also calculated the quantum amplitude with the complete normalization, 
up to the Kahler normalization. In this there is no contribution from the geometric 
distribution of the fields, since the contributions from the normalization and that from 
the dimensional reduction cancel. 

Besides the a'^/^ suppressions from the closed string couplings and oscillator normal- 
izations, the coefficient is given by ratios of products of gamma function, whose argument 
is of 0(0. 1), thus we expect a factor of 0(1) from the quantum amplitude. The dominant 
amplitude is classical one, which is exponentially suppressed. Thus it is easy to generate 
hierarchy of Yukawa couplings. However, from the top-down approach, it is not easy to 
locate the desired fields at the desired positions. 
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A Useful formulae 

The four point amplitude is described by the standard hypergeometric functions. The 
following relations 

F{a, b; c; z) =^^^^^^^^-^^F{a, h- a + h + I - c'l - z) 
r(c — a)T[c — b) 

+ - zr-'^~'F{c - a, c - 6; 1 + c - a - 6; 1 - z), (106) 

T[a)V{b) 

F{a, b; c; z) =(1 - 2)"-"-^F(c - a,c - b;c; z), (107) 
are useful. 

The classical action contains the following functions, 



Fo(l - X) = -(-1)-W^5(M, |)F(|, 1 - I + ^; 



X) 



P^(^x) = (^ — iy''2+k3)/N ^-l+{ki+k2)/N J^(^ki hz'^p(^hl. ]_ _ ^1 _|_ 



X] 



F,{1 -x) = (-l)-i+'^3/A^(i _ ^)-i+fe+A'3)/iV5(M, 1 _ M + I; 1 _ 

F3(a;) = 5(1, ^)F(f^, 1 - ^, M + (108) 
Fi'(x) =Fi(A;i^ Ar-A;i;a;), 

F!,{l-x)=F2{h^ N -h-l-x), 

where B{a,b) = r(a)r(6)/r(a + b) and it is the Euler beta function. For Fq and F3 
we have used above relations (11061) . (llOyp . as well as (|24|) . For ^ fcj = 2nN, the above 
functions become as 



Fo{l-x) 



sin(7rA;4/A^) 



sm(7rA;i/A') 

+ sin(vr(fe, + A:.)/iV) ^_.,^,,,^ _ ' 
sm(7rA;i/A') 



xi?(l-f,l-f)i?(M,M)-ii.3(^), (109) 

l-{k-i+ki)/N _ -^l-{k2+kz)IN -l-{k^+k2)/N 



Fl^il -x) = (_i)i-(fe3+'=4)/Af(i _ 



x5(l-f,l-f)i?(f,f)-iFo(l-x). 

For factorization of the amplitude, we need the asymptotic behaviors of hypergeometric 
functions in the z ^ 00 limit 

^/ , X ^,-.r(c)r(6-a) „ .,r(c)r(a-6) , , , 

F(a.6;c;l-.)^ ^i?iy-°i .-° + ^if';-^i ^-^ (111) 

r(o)r(c — a) r(a)r(c — 0) 
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For the higher order amphtude than four, we need a multivariable generahzation of 
hypergeometric function, called Lauricella D function [36]. It is defined as 



oo oo 



('^)miH \-mr(J^l)mi ■ ■ ■ {br)mr mi r. 

"3/1 • • • X 



(112) 

with \xi\ < 1 for all i. Here (a)„ is the Pochhammer symbol meaning 

r(a + n) 

(a)„ = — — = a{a + 1) . . . (a + n - 1), n > 0, 

la , , 

{a)n = 77 n<0. 

(1 - aj-n 

We can express the integration as [T7] 

N-l 

X F\^-^\h/N, 1 - fci/iV, . . . , 1 - fci-i/iV; fc./iV + h,+i/N- Xi,u . . . x,,L-i), 
pi =e-i-ik^/^+^)B{kN/N, ki/N) 

X Fif-')(A;^/iV, 1 - h/N, . . . , 1 - kN-2/N; ki/N + k^/N; 1 - Xa, . . . , 1 - XN-2), 

(114) 

where the conformally invariant cross ratio. 
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